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Abstract 

A generic class of metamaterials is introduced and is shown to exhibit frequency dependent 
double negative effective properties. We develop a rigorous method for calculating the frequency 
intervals where either double negative or double positive effective properties appear and show 
how these intervals imply the existence of propagating Bloch waves inside sub- wavelength struc- 
tures. The branches of the dispersion relation associated with Bloch modes are shown to be 
explicitly determined by the Dirichlet spectrum of the high dielectric phase and the generalized 
electrostatic spectra of the complement. 
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1 Introduction 

Metamaterials are new class of engineered materials that exhibit electromagnetic properties not 
readily found in nature. The novelty is that unconventional electromagnetic properties can be 
created by carefully chosen sub- wavelength configurations of conventional materials. The distinctive 
properties of metamaterials are derived from geometrically induced resonances localized to specific 
frequencies. These resonances are used to control propagating modes with wavelengths longer than 
the characteristic length scale of the material. Metamaterials are envisaged for several application 
areas ranging from telecommunication and solar energy harvesting to the electromagnetic cloaking 
of material objects. 

A generic metamaterial comes most often in the form of a crystal made from a periodic array of 
scatterers embedded within a host medium. The physical notions of frequency dependent effective 
magnetic permeability and dielectric permittivity are used to describe the behavior of propagating 
modes at wavelengths larger than the length scale of the metamaterial crystal. The past decade 
has witnessed the development and identification of new sub- wavelength geometries for novel meta- 
material properties. These include the simultaneous appearance of negative effective dielectric per- 
mittivity and magnetic permeability. Such "left handed media" are predicted to exhibit negative 
group velocity, inverse Doppler effect, and an inverted Snell's law |40j . The first metamaterial con- 
figurations imparted electromagnetic properties consistent with the appearance of a negative bulk 
dielectric constant [33]. Subsequently electromagnetic behavior associated with negative effective 
magnetic permeability at microwave frequencies were derived from periodic arrays of non-magnetic 
metallic split ring resonators |35| . Double negative or left handed metamaterials with simultaneous 
negative bulk permeability and permittivity at microwave frequencies have been verified for arrays of 
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metallic posts and split ring resonators [IS] . Subsequent work has delivered several new designs using 
different configurations of metallic resonators for double negative behavior [TSJ [221 Ull HSl [SOI [52] ■ 

Current state of the art metallic resonators do not perform well at optical frequencies and alter- 
nate strategies are contemplated employing the use of both metals and dielectric materials for optical 
frequencies |43| . For higher frequencies in the infrared and optical range new strategies for gener- 
ating double negative response rely on Mie resonances generated inside coated rods consisting of a 
high dielectric core coated with a dielectric exhibiting plasmonic or Drude type frequency response 
at optical frequencies |46], [47], [48]. A second strategy for generating double negative response 
employs dielectric resonances associated with small rods or particles made from dielectric materials 
with large permittivity, [27], [37], [41]. Alternate strategies for generating negative permeability at 
infrared and optical frequencies use special configurations of plasmonic nanoparticles [T] , [35] . The 
list of metamaterial systems is rapidly growing and comprehensive reviews of the subject can be 



Despite the large number of physically based strategies for generating unconventional properties 
the theory lacks mathematical frameworks that: 

1. Provide the explicit relationship connecting the leading order influence of double negative 
behavior to the existence of Bloch wave modes inside metamaterials. 

2. Provide a systematic identification of the underlying spectral problems related to the crystal 
geometry that control the location of stop bands and propagation bands for metamaterial 
crystals. 

In this article we provide such a framework for a generic class of metamaterial crystals made from 
non-magnetic constituents. The crystal is given by a periodic array of two aligned non-magnetic 
rods; one of which possesses a large frequency independent dielectric constant while the other is 
characterized by a frequency dependent dielectric response. In this treatment the frequency depen- 
dent dielectric response e_p is associated with plasmonic or Drude behavior at optical frequencies 
given bv [46 ] . [47 ] 



where lo is the frequency and is the plasma frequency [6] . Here we develop a rigorous method for 
calculating the frequency intervals where either double negative or double positive bulk properties 
appear and show how these intervals imply the existence of Bloch wave modes in the dynamic 
regime away from the quasi static limit, see Theorems 14.21 and 14.31 It is shown that these frequency 
intervals are explicitly determined by two distinct spectra. These are the Dirichlet spectrum of the 
Laplacian associated with the high dielectric rod and the electrostatic spectrum of a three phase 
high contrast medium obtained by sending the dielectric constant inside the high dielectric rod to oo. 
The electrostatic spectra is introduced in Theorem 12.11 and discussed in section [5] The methods are 
illustrated for ep given by (|l.ip , however they apply to dielectrics characterized by single oscillator 
or multiple oscillator models that include dissipation and are of the form 



where ujj are resonant frequencies, and 7j are damping factors. 

We start with a metamaterial crystal characterized by a period cell containing two parallel 
infinitely long cylindrical rods. The rods are parallel to the X3 axis and are periodically arranged 
within a square lattice over the transverse x = (a;i,a;2) plane. The period of the lattice is denoted 
by d. There is no constraint placed on the shape of the rod cross sections other than they have 
smooth boundaries and are simply connected. The objective is to characterize the branches of the 
dispersion relation for for H-polarized Bloch-waves inside the crystal. For this case the magnetic 



found in [H] and [33]. 



(1.1) 




(1.2) 
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field is aligned with the rods and the electric field lies in the transverse plane. The direction of 
propagation is described by the unit vector k ~ K2) and k = 27r/A is the wave number for a 
wave of length A and the fields are of the form 

^ iJ3(x)e'('^''^-''-*"/"), El = i;i(x)e'('='^-^-*"/"\ = £;2(x)e*('='^-^-*'^/'=) (1.3) 

where H^{x.), £'i(x), and £'2(x) are d-periodic for x in R^. Here c denotes the speed of light in free 
space. We denote the unit vector pointing along the X3 direction by 63, and the periodic dielectric 
permittivity and magnetic permeability are denoted by ad and /i respectively. The electric field 
component E = (i?i,i?2) of the wave is determined by 

ic 

E = 63 X VH3. 

Load 

The materials are assumed non-magnetic hence the magnetic permeability is set to unity inside 
the rods and host. The oscillating dielectric permittivity for the crystal is a d periodic function in 
the transverse plane and is described by aa — ad{'x./d) where arf(y) is the unit periodic dielectric 
function taking the values 

{Eh in the host material, 

ep = ep(w) in the frequency dependent "plasmonic" rod, (1.4) 
— £r/d'^ in the high dielectric rod. 

This choice of high dielectric constant ep follows that of [5] where €r has dimensions of area. 
Setting h'^{x.) = H^{-x)e^^^'^'^^ the Maxwell equations take the form of the Helmholtz equation given 
by 

-Vx- (a-i(^)Vx/i''(x)) =^/i^ inR2. (1.5) 

2 

The band structure is given by the Bloch eigenvalues ^ which is a subset of the parameter space 

-27r < fci < 2tt, -27r < ^2 < 27r}, with k = [kl + k^Y^^ and ki = kjk. This constitutes the 
first Brillouin Zone for this problem. 

We set X = dy for y inside the unit period Y = [0, 1]^, put /? — dkk and write u{y) = iJ3(dy). 
The dependent variable is written u'^{y) = h'^{dy) = ii(y) exp*^'^, and we recover the equivalent 
problem over the unit period cell given by 

-Vy(a^i(y)VX) = ^^' -^^Y. (1.6) 

To proceed we work with the dimensionless ratio p = d/y^, wave number t = i/e7fc and square 
frequency ^ = ^r^s-- The dimensionless parameter measuring the departure away from the quasi 
static regime is given by the ratio of period size to wavelength rj — dk — pr > 0. The regime i] > 
describes dynamic wave propagation while the infinite wavelength or quasi static limit is recovered for 
r] = 0. Metamaterials by definition are structured materials operating in the sub- wavelength regime 
< ?7 < 1 away from the quasistatic limit |35| . For these parameters the dielectric permittivity 

2/2 

takes the values ep — I — f!±i£_fL^ — — and is denoted by ap(y) for y in Y and ()1.6p is 

given by 

- Vy • {a;\y)Vyu'{y)) = p'^u''{y) in Y (1.7) 

The unit period cell for the generic metamaterial system is represented in Figure [TJ In what follows 
R represents the rod cross section containing high dielectric material, P the cross section containing 
the plasmonic material and H denotes the connected host material region. 

It is shown in Theorem 14.11 that the band structure for the metamaterial is characterized by a 
power series in 77 and is governed by two distinct types of spectra determined by the shape and con- 
figuration of the rods inside the period cells. The series delivers the explicit relationship connecting 
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Figure 1: Cross section of unit cell. 



the leading order influence of quasi static behavior, as mediated by effective magnetic permeability 
and dielectric permittivity, to the propagation of Bloch wave solutions inside metamaterial crystals 
made from sub- wavelength 1 > ry > structures, see Theorems 14.11 l42l and 14.31 

The relevant spectra for this problem is found to be given by the Dirichlet spectra for the 
Laplacian over the rod cross sections R together with a generalized electrostatic spectra associated 
with the infinite connected region exterior to the rods. These spectra provide two distinct criteria 
that taken together are sufficient for the existence of power series solutions see, Theorem 14.21 In 
what follows the power series is developed in terms of a hierarchy of boundary value problems posed 
separately over the domain R and the domain exterior to the high dielectric rod Y\R. The existence 
of solutions for the boundary value problems inside the high dielectric rod R is controlled by the 
Dirichlet spectra see, section [3l Existence of solutions for boundary value problems exterior to R are 
determined by the generalized electrostatic spectrum see. Theorems 12.11 and 12.21 The electrostatic 
spectra is identified and is shown to be given by the eigenvalues of a compact operator acting on an 
appropriate Sobolev space of periodic functions see, the discussion in section [5] and Theorem 15.21 

The generic class of double negative metamaterials introduced here appears to be new. The 
motivation behind their construction draws from earlier investigations. The influence of electrostatic 
resonances on the effective dielectric tensor associated with crystals made from a single frequency 
dependent dielectric inclusion is developed in the seminal work of [5], see also the more recent 
work |38| in the context of matematerials. These resonances are responsible for negative effective 
dielectric permittivity. In this context we point out that the sub-wavelength geometry introduced 
here is a three phase medium and the categorization of all possible electrostatic resonances requires a 
different approach see section [S] On the other hand Dirichlet resonances generate negative effective 
permeability inside high contrast non-dispersive dielectric inclusions, this phenomena is discovered in 
[IB] , [5] , [TU] , see also the mathematically related investigations of [HI] , [H] , [12] and [IS] . Motivated 
by these observations we have constructed a hybrid composite crystal that combines both high 
dielectric inclusions and frequency dependent inclusions for generating double negative response 
from non-magnetic materials. 

The power series approach to sub-wavelength 77 < 1 analysis has been developed in [23] for 
characterizing the dynamic dispersion relations for Bloch waves inside plasmonic crystals. It has 
also been applied to assess the influence of effective negative permeability on the propagation of 
Bloch waves inside high contrast dielectrics [24], the generation of negative permeability inside 
metallic - dielectric resonators ;39 , and for concentric coated cylinder assemblages generating a 
double negative media [TT] . 

We conclude noting that earlier related work introduces the use of high contrast structures for 
opening band gaps in photonic crystals, this is developed in [T^, [5^, and [H]. For two phase 
high contrast media integral equation methods are applied to recover dispersion relations about 
frequencies corresponding to Dirichlet eigenvalues [3] and [4]. The connection between high con- 
trast interfaces and negative effective magnetic permeability for time harmonic waves is made in 
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[30) . More recently two-scale homogenization theory has been developed for three dimensional split 
ring structures that deliver negative effective magnetic permeability [S], |13) . For periodic arrays 
made from metal fibers a homogenization theory delivering negative effective dielectric constant 
[7] is established. A novel method for creating metamaterials with prescribed effective dielectric 
permittivity and effective magnetic permeability at a fixed frequency is developed in |32) . 



2 Background, basic theory, and generalized electrostatic 
resonances 

In this section we outline the steps in the power series development of Bloch waves and establish 
the mathematical foundations for establishing existence of power series solutions. We introduce the 
function space H^^^iY) defined to be all F-periodic, complex valued square integrable functions 
with square integrable derivatives with the usual inner product and norm given by 

(Vu • VtJ + uw) and = (m, u)y . (2-1) 

We also introduce the space H^^^iY \ R) given by all y-periodic, complex valued square integrable 
functions with square integrable derivatives and the inner product and seminorm 



{u,v)^i Vw • Vwdy and ||w|| (2.2) 
Jy\r 



IY\R 

The variational form of (jl.7p is given by 

a-^Vu'^-Vi= [ ^u'^v (2.3) 
/y Jy c 

for any v — v{y)e^'^''^P'^ , where v £ H^^^iY) . On writing m'' ~ u{y)e^'^''^Py and setting rj — rp we 
transform (|2.3p into 



T^(^ — er^)(V + irjk)u ■ (V + irjHi)v + I t^^(V + ii]k)u ■ (V + i'qk)v 

c Jp 

- er^)(V + ir]k)u ■ (V + iT]k)v = / rf-i^i^ er^)uv (2.4) 
j ji c Jy c 

We introduce the power series 



oo 

U — 

m=0 



V'^Um (2.5) 

oo 



m=0 



where Um belongs to Hp^^{Y). In view of the algebra it is convenient to write Um — i'"uoV'm where 
Wg is an arbitrary constant factor. 

We now describe the underlying variational structure associated with the power series solution. 

Set 



z = ep\^o)= 1- 



Co ; 

and for u, v belonging to Hp^^.{Y) we introduce the sesquilinear form 

B:,{u,v) ^ I Vu-ViJdy+ / zVu-Vwdy. (2.7) 
Jh Jp 
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Here Y\R = HU P and the form Bz{u, v) is well defined for functions in Hp^^{Y \ R). Substitution 
of the series into the system (j2.4p and equating like powers of 77, produces the infinite set of coupled 
equations for m = 0, 1, 2 . . . given by 

m—1 m — 1 



n III ± IIL J. 

^o^^p\^oy / [ V(-^)'6VV'm-^•V^; + K• V(-i)'6(V'm-i-iVw-V?/'™_i_iTJ) 

■^y\R' 1=1 1=0 

m-2 ^2 ^ 

(-i)'6V'm-2-il'] - io^<^p^{io)T'^<^r^ / • (V-m-lVU- Vl/jm-iv) " ■0m-2t^] 

1=0 ^ 



m — 3 



'(Co) / [ i-^yWrn-2-l • ViJ + K ^ ("^ )'6 (V-m-S-i ViJ - V7/>™_3-it^) 

•^■'^ ;=o 1=0 

™-4 „ ^2 

y2{-iy£.l1p7n-4-lv] +Ccr^ep^(Co) / er^[V-0m-2-VlJ+K('0m-3VlJ-V?/'m-3^') + '/'m-4w] 

„ 771—2 ^ 2 
C(7'S '(Co) /[EE Cm-2-ienV'i-n»'""""^+ e,-^ E ("O'^ V'm-2-i^7] 



i=0 n=0 i=0 



= 0, for ah w in iJpi,^(r). (2.8) 

Here the convention is ipm — for to < 0. 

The determination of {i^m}m=o proceeds iteratively. We start by determining ■0o on y \ R, this 
function is used as boundary data to determine ipo in R from which we determine tpi on Y \ R 
and the full sequence is determined on iterating this cycle. The elements are recovered from 
solvability conditions obtained by setting u = 1 in (|2.8p and proceeding iteratively. The complete 
algorithm together with explicit boundary value problems necessary for the determination of the 
sequences {ipm}m=o^ {^m}m=o described in section [3] and Theorem 13. II 

The existence theory for the solution of the sequence of boundary value problems is based on 
the sesquilinear form Bz{u,v) defined for functions u and v belonging to Hp^^{Y \ R)/C Here 
Hp^^{Y\R)/C is the subspace of functions u belonging to Hp^^{Y\R) with zero mean Jy^^j^ udy = 0. 
This space is a Hilbert space with inner product (12. 2p . Although in this treatment ep is given by 
()l.ip we are motivated by the general case (|1.2p and proceed in full generality allowing for the 
possibility that z can lie anywhere on the complex plane C including the negative real axis. Thus 
for each z in C we are required to characterize the range of the map u i-> Bz{u, •) viewed as a linear 
transformation mapping u into the space of bounded skew linear functionals on H^^^{Y \ R)/C. 
This is linked to the following eigenvalue problem characterizing all pairs XinC, ip in Hpg^(Y\R)/C 
that solve 

-i J^y^.yydy+l J^\/ij.\/vdy={Xi;,v), (2.9) 

for every v in H^^^{Y \ R)/C. Inspection shows that for z = {X + l/2)/(A - 1/2) that B^iip, v) = 0, 
for every v in Hp^j,{Y \ R)/C. In other words the kernel of the operator Tz is nonempty for z ~ 
(A + l/2)/(A — 1/2). The eigenvalues A will be referred to as generalized electrostatic resonances. 
In what follows we show that Tz is a one to one and onto map from H^^^ (Y \ R)/C into the dual 
space provided that z {X + l/2)/(A - 1/2). 

The generalized electrostatic resonances are characterized by introducing a suitable orthogonal 
decomposition of H^^^iY \ R)/C We introduce the subspace Hq{P) given by the closure in the 
norm of smooth functions v with compact support on P and the subspace Hq p^^(H) given by the 
closure in the norm of all periodic continuously differentiable functions with support outside P. 
Extending elements oi Hq p^j.{H) by zero to F \i? delivers Wi C Hp^j.(Y\R), extending elements of 
H^iP) by zero to F \ i? delvers W2 C H^^^{Y \ R). Define W3 to be all functions w in H^^^iY \ R) 
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for which the boundary integral J^p w dS vanishes and that belong to the orthogonal complement 
of Wi U W2 with respect to the inner product (|2.2p . Its easily verified that Wi, W2, and are 
pairwise orthogonal with respect to the inner product (j2.2l) and 

H^^^{Y\R)/C^Wi®W2®W3®C, (2.10) 

where the constant part of a function u belonging to this space is uniquely determined by the 
condition Jy^n udy — 0. 

The following theorem describing all eigenvalue eigenfunction pairs is established in section [5l 

Theorem 2.1. The eigenvalues for ()2.9p are real and constitute a denumerable set contained inside 
[—1/2,1/2] with the only accumulation point being zero. The eigenspaces associated with A = 1/2 
and A = —1/2 are Wi and W2 respectively. Eigenspaces associated with distinct eigenvalues in 
(—1/2, 1/2) are finite dimensional, pairwise orthogonal, and their union spans the subspace W3. 

We denote the denumerable set of eigenvalues for (|2.9p by the sequence {Xi}°Zi. The orthogonal 
projections onto Wi and W2 are denoted by Vi and 1^2- The orthogonal projections on to the finite 
dimensional eigenspaces associated with the eigenvalues Ai in (—1/2, 1/2) are denoted by Vx-. For 
z {\i + l/2)/(Ai — 1/2) the following existence theorem holds. 

Theorem 2.2. Suppose z ^ {Xi + l/2)/(A, - 1/2) for X, G [-5, 5] then 

• For any F G [Hp^j,(Y \ R) /C] * such that F(v) — for constant v, there exists a unique solution 
u e Hp^^(Y \R)/C of the variational problem Bz{u,v) ~ F{v) for all v G H^^^{Y \R)/C 

• The transformation T^ from Hp^^{Y\R)/C onto itself has the representation formula given by 

T,=Vl+zV2+ J2 + (2-11) 

-i<A„<i 

with inverse 

T-' =Vi+z-'V2+ J2 (1 + (^-1)(^-A„))-^Pa„, (2.12) 

-i<A„<i 

and Bz{u, v) = {TzU, v) for all u, v in H^^^{Y \ R)/<C. 

This theorem is proved in section |6l 

In the following sections we present the sequence of boundary value problems for determining ■0™ 
in y \ i? and R together with the sequence of solvability conditions characterizing In this section 
we use the complete orthonormal systems of eigenfunctions associated with electrostatic resonances 
and Dirichlet eigenvalues to explicitly solve for fields ipo in Y and ipi uiY\R and provide an explicit 
formula for ^o- In the following section the boundary value problems used to determine ipm in y \ i?, 
m > 2 and ipm hi R for m > 1 together with solvability conditions for for m > 1 are shown to 
be a well posed infinite system of equations see Theorem 13.11 

Applying the convention = for m < in (|2.8p shows that ipQ is the solution of 

Bz{^po,v)^0, forallwiniJi(r\i?). (2.13) 
Prom Theorem 12.11 and ()2.10|) we have the dichotomy: 

1. ^0 satisfies ep^(^o) = (Xi + l/2)/(Ai — 1/2) and tpo is an eigenfunction for (|2.9p . 

2. ^0 satisfies ep^(^o) i^i + l/2)/(Ai — 1/2), i = 1,2, . . . and tAo — constant. 
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In this article we assume the second alternative. Subsequent work will investigate the case when 
the first alternative is applied. For future reference the condition ep^(^o) 7^ {K + l/2)/(Ai — 1/2) is 
equivalent to 



e 

Co^C:, G = (A. + i)^ ,0<C.<^* = 1,2,.... (2.14) 
Restricting to test functions v with support in R in (|2.8p we get 

(VV'o • Vw - Co^ow) rfy = 0. (2.15) 



From continuity we have the boundary condition for ipQ on R given by ipQ = const. We denote the 
Dirichlet eigenvalues for i? by ia, , j = 1, 2, . . .. Here we have the alternative: 

1. If ^0 is a Dirichlet eigenvalue i/i of —A in R then tpoly) = for y in F \ i?. 

2. li ^0 i^i, i — 1,2, .. . then V'o is the unique solution of the Helmholtz equation (|2.15p and 
tjjo = const, in y \ i?. 

In this treatment we will choose the second alternative 7^ ^'i- The case when the first alternative 
is chosen will be taken up in future investigation. Since uq = UqiPo where Ug is an arbitrary constant 
we can without loss of generality make the choice ipo = 1 for y in 1" \ i?. Since 7^ ^ii i = ^, ■ ■ ■ 
and ipo = 1 in Y \ R a. straight forward calculation gives 4'o in R in terms of the complete set of 
Dirichlet eigenf unctions and eigenvalues: 

00 

^0 = > , z m R. (2.16) 



71=1 



Note here that /i„ denote the Dirichlet eigenvalues of — A in i? whose eigenfunctions 0„ have nonzero 
mean, < >_r= (j)n{y)dy 7^ 0. The Dirichlet eigenvalues associated with zero mean eigenfunc- 
tions are denoted by and = {Mn}^=i U Wn}n=i- 

To find tpi in F \ R, we appeal to (|2.8I) with i/jq = 1 in Y \ R to discover 



B,ii;i,v) = - [ k-Vv- j ep\S,o)k-Vv V v e hI^,{Y) (2.17) 

J H Jp 



It follows from Theorem 12.21 that the problem has a unique solution subject to the mean-zero 
condition: ^y\r. V'l = provided that ^0 7^ Cii * = 1: 2, . . .. We apply the decomposition of H^^^ {Y \ 
R) /C given by (j2.10p and represent ipi in terms of the complete set of orthonormal eigenfunctions 
{?Aa„} C W3 associated with —1/2 < A„ < 1/2 together with the complete orthonormal sets 
of functions for Wi and W2, denoted by {V'nlS^Li and {i^n}^=i respectively. A straight forward 
calculation gives the representation for -0i in F \ i? 

l + (.p\Co)-l)(l-A„)J^^"-^g"^-^"' ^'-''^ 




with 



■ / ^"^^^^ '^y^ "a„ = ^ ■ /p ^V'a„ c^y, and = k- Vip^ dy. (2.19) 

JH " JH 

Setting V — 1 and m = 2 in (j2.8l) we recover the solvability condition given by 

[ [~k- CoVVi + Co] - rhr'^ [ {-k\/^lj, + 1) (2.20) 

(CoV'O - Er^CoV'o) 
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Substitution of the spectral representations for ipi and i/jq given by (|2.16l) and (I2.18P into (I2.20p 
delivers the quasistatic dispersion relation 



Co 



2 -2 

T n 



'eff 



(Co) 



(2.21) 



where the effective index of diffraction 
written 



depends upon the direction of propagation k and is 



(2.22) 



The frequency dependent effective magnetic permeability fJ.eff and effective dielectric permittivity 
Ce/ / are given by 



Me// (Co) 



-00 



Co 



(2.23) 



and 



^e//(Co)K-K 



(/-P)k- Kdy + 



Co 



H 



to ;3- 



E 

-l/2<A,,<l/2 



[^0 - 



C At A>, 



Co- 



Co - S/i 



V 



(?.24) 



where and 9p are the areas occupied by regions H and P respectively. The first term on 
the right hand side of (|2.24l) is positive and frequency independent, ft is written in terms of the 
spectral projection P of square integrable vector fields over H onto the subspace of gradients of 
potentials V in Hqp^^^H). Here Jj^Pk ■ kdy < 9h and Per = Y^^=i (/f/^V'n ■o'rfy) VV'i with 
l"^i,nP = Ih^'^ ' ^dy. The poles Sh are the subset of values Ch for which at least one of the 
weights o'^^ and a'^^ are non zero. The values of C^h for which both weights vanish are denoted 
by sJi and {Chj^Li — {s/i}/^Li ^ The graphs of ^leff and ^^ffk ■ k as functions of are 

displayed in Figure [5] and Figure [31 Here the intervals a' < < b' and a" < Co < b" are the same 
in all graphs. 

Me// 



a' b' /a" h"' 



Ml 



Hi 



Co 



M3i 



Figure 2: The relation between /ie// and 



For future reference is is convenient to write the dispersion relation explicitly in terms of ^^eff 
and Ce// and 



Me//(Co)Co = T'^eJf{^o)K ■ k. 



(2.25) 



Figure 3: The relation between ^eff^ ' ^^'^ ^o- 



The branches of the quasistatic dispersion relation (r, k) i— > a-re controlled by the poles and zeros 
of fieff and e~jj: explicitly determined through the Dirichlet spectra and generalized electrostatic 
resonances. 



3 Solution of higher order problems 

Having determined ipo in Y, i/ji in y \ R, and we now provide the algorithm for determining the 
rest of the elements in the sequences {ipm}m=0' {^m}m=o- The algorithm is summarized in Theorem 

Step I. Solution of ipm in R for m > 1 . 

We restrict the trial space to test functions v with support in R and put m i— > to + 2 in (|2.8p . 
We decompose ipm according to 

V-m = V'm + {-iTim^* (3.1) 

and substitute into (|2.8p . This decomposition is chosen such that V'm depends on ^„ and ipn inside 
i? for n < TO — 1. The function depends only upon and ijjQ in R. The function '^rn solves the 
Dirichlet boundary value problem with V'r?i|a/?,- ~ ''Pm\dR+ a-nd 

j i^^rn ■ VtJ - ^o^mxj) dy ^ j F- Vvdy + J Gvdy, (3.2) 



where 



m— 1 m — 1 



F = -eo"'s'(^o)[E(-*)'^'^^™-' + ^(E(-*)'^'^™-i-')] 

1 = 1 1=0 
^2 

Ho'^p\^o)er^k,Pm-i, (3.3) 

and 

m — 1 m — 2 

G = ^O^'S '(^0)(Ai • (-*)'aVVm-l-i + E (-*)'6^m-2-0 

1=0 1=0 

-io^e;\(o){k-^4'm-l+^m^2) (3.4) 
rn— 1 I 2 TTi — 1 m — 1 

1=1 n=0 1=1 n=l 

This boundary value problem has a unique solution for 7^ ^'i: * = 1, • . •. The function ip^, is the 
solution of 
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\'(p*\dR- = 

The explicit representation for ?/'* is obtained using the Dirichlet eigenfunctions on R and is given 
by, 

oo , 

Step II. Solution of iprn inY\R for m>2. 
We decompose V'm 

V'm = V4 + {-iT-^U^i^ (3.7) 

and substitute into (|2.8p . This decomposition is chosen such that depends on ^„ for < ri < 
m — 2, the functions '(/'„,n<m— liny\i? and the functions ij^m n < m — 2 m. R. The function if) 
depends only on and ipi in F \ R. 

For m > 2, ?/;^ vaY \ R subject to the mean zero condition ^y\r c^y = is the solution of 

r'^B^i^P'^, v)^ ( (Fi • VtJ + Giv) dy + [ (F2 • Vv + G2V) dy, for all v in iJpV(r), (3.8) 
Jy\r Jr 

where 

m— 2 m— 2 

^1 = Co"iepi(Co)r2[EH)'6VVm-i + (E(-*)'^'V'™-i-0«] 
1=1 1=0 
^2 

-fo'^S ^(fo)T^er-^XH«^V'm-i, (3.9) 



Gi = ^o'^p'i^oy [f^ ■ (-*)'6(-vv^™-i-i) - 

1=0 1=0 

-^o^'^P^i^o)r^er-^XH[K ■ (-VV'm-i) - ipni-2] 

m — 2 I 2 m — 2 

+eo"'s '(^O) [ E E e™-2-ienV'/-n*'-"-" " ^ ^ , (3.10) 

in — 2 m ~ 3 



-^(7'^ '(Co) [ E (-^)'6vv™-2-i + «(E (-o'av'm-s-i)] 

^2 

+eo"'S '(^o)er.^ [VVm-2 + Ac(Vm-3)] , (3.11) 



and 



m — 3 



G2 = -Co"^6;l(eo)[ + ^(E(-*)'6(-VV'm-3-/)) 

/=0 

m-4 ^2 

- E (-«)'C;V'm-4-(] + C(7^ep ^(Co)er^ [k ■ (-VV'm-s) + V-m-J 
(=0 

m—2 i 2 m — 2 



i=0 n=0 
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where xh is the indicator function of the host domain H taking 1 inside and zero outside. Here we 
have used the identity -00 (y) — 1 > for y Y \ R in determining the formulas for the right hand side 
of (03). 

Setting w = 1 in p.8p dehvers the solvabihty condition 



Gidy+ / G2dy = 0. (3.13) 

Y\R Jr. 

The condition 

- V • F2 + G2 = for y in i? (3.14) 



follows from the solution of Step I see, (13. ip . (|3.2p and (|3.5p . Integration by parts on the left hand 
side of p.Sp together with p.l4p transforms (|3.8p into the equivalent Neumann boundary value 
problem for ip'^, m > 2, given by 

T^B^{iP'^,v)^ [ {Fi ■ S/v + Giv) dy + [ F2 ■ nv ds, for all v in H^^^{Y \ R), (3.15) 

Jy\R JdR. 

where n is the outward directed unit normal to dR. Finally we observe that the solvability condition 
for the Neumann problem p.lSp given by 



Gidy + / F2-nds = 0, (3.16) 

Y\R JdR 

follows immediately from (|3.13p and p.l4p . This Neumann problem satisfies the hypotheses of 

Theorem 12.21 and we assert the existence of a solution for p.isp , uniquely determined by the 

2 

condition /^^^ dy = 0, provided that ^0 7^ Ci, < d < i = 1, 2, . . .. 

The field ip solves 

(3.17) 

for all trials v G H^^^ {¥ \ R) . The solution i/j is represented explicitly in terms of the eigenvectors 
associated with the generalized electrostatic resonances. A straight forward calculation gives 



l/2<A„<l/2 Co — (1/2 + Xn)-^ 



Co - Si 

\ J 



(3.18) 



Step III. Solution of S,m for m > 1. 

We apply the decomposition p.ip to ipm-2 in R and (|3.7p to i/im-i in F \ i?. Substitution of the 
decompositions into the solvability condition p.l3p and applying the explicit spectral representations 



12 



(|3.6p . (j3.18p for -0* and i/j we obtain the recursion relation for determining ^i, ^2, • ■ • given by 

771 — 3 771-3 



1=1 1=0 

^« /=0 1=0 

„ 771 — 3 ^ 771 — 3 p p 

-/[EE ^'~""'"e™-2-ienV';-n - E (-*)'eo6V'™-2-;] - / eo^m-2 - / Co' 



2 ; 
V'm-2 



I Er— > (-«) 6'/'m-2-i + / ^r^^O'4>rn-2 + I 6^— ^oV'm-2- (3.19) 

Jy c j_fj., c Jy\r c 



Noting that (I3.ip and p.7p can also be applied to the lower order terms on the right hand side ol 
(|3.19p we see that (I3.19|) determines ^m-2 from ^„ in F \ i?, and ijjn in R with < n < m — 3. 
Here the "solvability matrix," Q depends explicitly on and is given by 

Q{io) = OH + ep 

« oo 

+T^[ / P^- «dy- 5^(5^(eo)|aL +^p'(?oXn + 25„(eo)« +<)« +ep'(eo)aL)] 

-^H n=l 

2 2 oo / / \2 

-2(^0 + ^)/.e//(eo) + eo(eo ~ ^) E r"^^ (3-20) 

2 2 

where gn{Co) = (^o - ^^)/(^o - Sn) and Cp^^o) = ^o/(Co - ^^)- The set of poles and zeros for 
G{^o) are explicitly controlled by the generalized electrostatic resonances and Dirichlet spectra. The 
zeros of G are denoted by 7„, n = 1, . . .. It is now evident that the recursion relation p.l9p holds 
provided that ^n, Co 7^ Cn, and <^o 7^ 7n- 

Denote the union of the Dirichlet spectra, the electrostatic resonances, and zeros of Q by 

U ^ U U {7,}°li (3.21) 

where C M+, {lj}f=i C M+, and {Cj}^i C [0,^]. We collect results and state the 

following theorem. 

Theorem 3.1. //^q belongs toR+\L/ then the sequences {V'm}m=o ^ Hp^^(Y), {S,m.}m=o <^ *C ^^list 
and are uniquely determined. 

Proof. The functions ?Ao in i?^(F) and -01 in iJ'^CK \ R) with jy^^ipidy — exist as does Co- 

Application of Step I uniquely determines ^"1 in R. Application of Step II uniquely determines ?/;2 
in y \ i?. Step III uniquely determines Ci and we apply p.ip and p.7p to recover -01 in if^(y) and 
02 in H^{Y \ R) with /^^^ -4,2 dy = 0. 

We now adopt the induction hypotheses H„, n > 2: 
There exist 

1) 0„ e H^,,(Y \ R) with Jy^^ V« = ; 

2) 0j e H^eri^) for < j < n - 1 with /^^^ = 0, 1 < j < - 1; 

3) Cn-i e C. 

Application of Step I uniquely determines 0„ in R. Application of Step II uniquely determines tpn+i 
in y \ i?. Step III uniquely determines Cn and we apply p.ip and p.7p to recover 0„ in H^{Y) and 
Tpn+i in i7i(F \ i?) with /^^^ Vn+i dy = 0. □ 
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4 Band structure for the metamaterial crystal 



In this section we present explicit formulas describing the band structure for the metamaterial 
crystal in the dynamic, sub-wavelength regime 1 > t] > 0. The formulas show that the frequency 
intervals associated with pass bands and stop bands are governed by the poles and zeros of the 
effective magnetic permittivity and dielectric permittivity tensors. The poles and zeros are explicitly 
determined by the Dirichlet spectra of R and the generalized electrostatic resonances Y \ R, see 
(|2.23p and (|2.24p . For the general class of parallel rod configurations treated here the pass bands 
and stop bands can exhibit anisotropy, i.e., dependence on the direction of propagation described 
by k. The anisotropy is governed by the projection of k onto the eigenfunctions associated with the 
generalized electrostatic resonances given by (|2.19p . In what follows pass bands are explicitly linked 
to frequency intervals and propagation directions for which both the effective magnetic permittivity 
and dielectric permeability have the same sign. This includes frequency intervals where effective 
tensors are either simultaneously positive or negative. 



r 




Figure 4: The leading order dispersion relation over two selected intervals /„/ and /„" . 



We begin by identifying the locations of the branches for the dispersion relation associated with 
the metamaterial crystal. We introduce the union of open intervals 1J„ 0„ on the positive real axis 
M+ obtained by removing the points {Cjjf^i, {i^jjfLi, {ijj'jLi, {s*j}j^i and Here s* is the 

zero of t~j'j{^a)k ■ k between Sj and Sj+i and /i* is the zero of fj,effi£,o) between fij and fJ-j+i- The 
leading order dispersion relation (j2.25p implicitly defines the map (r, k) i— > ■ We denote the branch 
associated with 0„ by = ^q"''(t, k) . Let /„ be an open interval strictly contained inside On- For 
this choice /„ C 0„ does not intersect the union of the Dirichlet spectra, generalized electrostatic 
spectra and the zeros of the solvability matrix Q, see (|3.21l) . 

The band structure for the metamaterial in the dynamic sub-wavelength regime 1 > ry > is 
given by the following theorem. 

Theorem 4.1. The dispersion relation for the metamaterial crystal contains an infinite sequence 
of branches with leading order behavior given by the functions = ^q(t, k). For r such that 
belongs to /„ there exists a constant R depending only on In such that for < p < R, the branch of 
the dispersion relation for the metamaterial crystal is given by 

oo 

^ = ^l}{r,k) + Y,{pry^r, (4.1) 

1=1 
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for 

{-27r < Tpki < 27r, -27r < rp^a < Stt}. (4.2) 

Here the higher order terms ^" are real and are uniquely determined by according to Theorem 

roi 

The power series representation for the transverse magnetic Bfoch wave (|1.3p is given by the 
foUowing theorem. 

Theorem 4.2. For In and R as in Theorem \4-l\ and for t such that belongs to In there exist 
transverse magnetic Bloch waves given by the expansion 

Hs^^Uq ^Mx/d) +Yy.PTyi^Mx/d)^ exp^i (kk- x-t'^'j^ (4.3) 

where the series 

oo 

My) + Y.'^pryi'My) (4.4) 

1=1 



is summable in H^(Y) for < p < R and ^ — 




Theorem 14.11 expUcitly shows that the leading order behavior determines the existence of pass 
bands or stop bands when p = dj ^JT^ is sufficiently small. With this in mind we can appeal to (I2.25|) 
and state the following theorem. 

Theorem 4.3. For p > sufficiently small: 

• There exist propagating Bloch wave solutions along directions k over intervals In for which 
PeffiCo ) o,''T'd ^effi^o )^ ' ^fl'f'c the same sign. 

• Intervals In for which PeffiCo) '"^'^ ^eff i^o)'^ ' have the opposite sign lie within stop bands. 

Leading order behavior for pass bands associated with double negative and double positive be- 
havior are illustrated in Figure 2] for two intervals /„' and Here the effective properties are both 
negative over the interval /„/ while they are both positive over /„" . 



5 Generalized electrostatic spectra 

In this section we establish Theorem l2.1l and characterize all pairs (A, ■(/;) in C x II^^^{Y\R) satisfying 
(|2.9p . We start by recalling the bilinear form in (12.91) and forming the quotient 

(u,m) 

From (|5.ip it is evident that A lies inside [—1/2,1/2]. It is easily seen that the solutions (A, ^/i) of 
(12. 9p for the choices A = ^ and A = — ^ correspond to ifi in Wi and W2 respectively. Only the 
element of W3 satisfies (|2.9p for the choices A = i and A = — i. 

We now investigate solutions (A, iji) for ij) belonging to W^^. The bilinear form (|2.9p defines a map 
T from W3 into the space of skew linear functional on Here W3 is a Hilbert space and 



(Tu,v) = — - j Vu-Vvdy+— j Wu-Wvdy, 
"2 Jp 2 Jfj 



(5.2) 
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for u and v in W3. In what follows we show that T is a compact map from W3 onto W3 with 
eigenvalues contained inside the open interval (—1/2,1/2). We do this by providing an explicit 
representation for T given in terms of a composition of single and double layer potentials. 

Set B = \Jnf=z,^{{Y \ R) + n) and introduce the Green's function G'(x,y) on B such that: G is 
separately periodic in x and y with period Y, G is in each of the variables x and y for x 7^ y, 
and 

AxG(x, y) = ^ ,5(x - y + n) - 1 in B, a„^G(x, y) = jS, for x e U„ez2 {dR + n), (5.3) 

where (5x is the Dirac delta function located at x and \R\ and \dR\ are the area and arc-length of R 
and dR respectively. We can find G(x, y) as a sum of the periodic free space Green's function F for 
the Laplacian and a corrector (jf: G(x, y) = F(x, y) + (/()*(x,y) , where 

nGZ2\{0} ' ' 

(/)*(x, y) is periodic in x and y with period F, G^ in x, y, and solves 
Ax</)*(x,y) =0 forx, yinB, 

dnA*{^^y)\+ = -dn^F{^iy)\+ + -p^jfor x on (9i?, 
aji \aR \dR\ 



and / (/)*(x,y)d5x = 0. (5.5) 

JdR 

Here the subscript indicates traces of functions on dR taken from the outside of R. 

The space of mean zero square integrable functions defined on dP is denoted by {dP) and for 
(j) in L'^{dP) we introduce the single layer potential Sp given by 

~Sp{cj,) = f G(x,y)0(y)d5y, (5.6) 

JdP 

and the modified single layer potential Sp mapping L^{dP) into given by 



Sp{4>) = Sp{4')-\dP\-^ Sp{4>){y)dSy. (5.7) 

JdP 

On restricting S'p(0)(x) to x on 9P, one readily verifies that Sp is self-adjoint with respect to the 
[dP) inner product and is a bounded one to one linear map from Ll{dP) onto H^^^{dP)nLl{dP). 
The inverse denoted by Sgp : H^/^{dP)r\L^{dP) — ?► L^{dP) is linear and continuous. We introduce 
the trace operator 7 mapping W3 into H^/^{dP) n LQ{dP), this map is bounded and onto 
Collecting results we have the following 

Lemma 5.1. Sp is a one to one, bounded linear transformation from L^{dP) onto W3 with Sp^ : 
W3 — > Lq^OP) linear and continuous given by Sp^ — Sgp'f. 

Proof. Given u in W3 consider its trace on dP given by 7W — u^^^ and u^^^ belongs to H^/^{dP) n 
L'^{dP). For X in y \ i? set w(x) — S'p(5gp(u|3p))(x). Since w belongs to W3 the difference w — u 
also belongs to W3 . Noting further that the traces of w and u agree on dP we conclude that w — u 
also belongs to Wi ® W^2- However since W3 = {Wi ® VF2)"'" it is evident that w — u = and we 
conclude that Sp^ = Sgp^. □ 

The outward pointing normal derivative of a quantity q on 9P is denoted by and the jump 
relations satisfied by the single layer potential are given by 

Sp{^):+ ^Sp{4,).- , on dP,^^f^\ ^^±h + K*p{cp), ondP (5.8) 
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where the double layer Kp is a bounded linear operator mapping Lq^OP) into LQ{dP) defined by 

Km = I ^^^P^mdSy, (5.9) 
JdP onx 

and the subscripts + and — indicate traces from the outside and inside of P respectively. Here 
is a continuous kernel of order zero and it follows ^22) that K*p is a compact operator mapping 
Ll{dP) into LQ{dP). Now we identify the transform T defined by ()5.2p . 



Theorem 5.2. The linear map T : W3 — >■ 14^3 defined by the sesquilinear form (|5.2p is given by 

T^SpK*pSp\ (5.10) 

and is a compact bounded self-adjoint operator on W3 with eigenvalues lying inside (—1/2, 1/2). 

Since T is self-adjoint, bounded, and compact: the spectrum is discrete with only zero as an accu- 
mulation point, the eigenspaces associated with distinct eignevalues are pairwise orthogonal, finite 
dimensional, and their union is W3. Theorem 12.11 follows immediately from Theorem 15 .21 noting that 
the choices A = 1/2, A = — 1/2 in (j2.9|) correspond to (j) belonging to Wi and W2 respectively. 
We now prove Theorem 15.21 



Proof Step I. 

We establish (|5.10p . For u, v, consider 

{SpK;Sp\u),v)^ [ VSpK;Sp\u)-Vvdx. (5.11) 
Jy\r 

Integration by parts gives 

' d{SpK;Sp\u)i 



{SpK;Sp\u),v) 



dP 



vdS, (5.12) 



+ 

where [q]^ = q\gp- — q\gp+- Applying the jump condition (|5.8p we get 



{SpK;SpHu),v) = - K;Sp\u)vdS, (5.13) 

JdP 



applying (15.81) again gives 



_i I dSpSp^u) ^ , I dSpSp'ju) ^ 



1 du 1 du 



(5.14) 



Step I now follows on substitution of (j5.14l) into (j5.13l) and integration by parts. 

Step II. The remaining properties of the operator T follow directly from the representation T = 

SpKpSp^ and the properties of Sp and Kp. □ 

We conclude noting that the spectrum of T is the same as the spectrum of Kp. This is stated 
in the following Lemma. 

Lemma 5.3. Tu — Xu if and only if X corresponds to an eigenvalue of Kp. 

Proof. If a pair (A,u) belonging to (—1/2,1/2) x W3 satisfies Tu = Xu then SpKpSp^u = Xu. 
Multiplication of both sides by Sp^ shows that Sp^u is an eigenfunction for Kp associated with 
A. Suppose the pair {X,w) belongs to (—1/2,1/2) x L^^dP) and satisfies KpW = Xw. Since the 
trace map from W3 to H^^^{dP) O L^{dP) is onto then there is a m in W3 for which w ~ Sp^u and 
KpSp^u = XSp^u. Multiplication of this identity by Sp shows that u is an eigenfunction for T 
associated with A. □ 
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6 Existence for exterior problems with a dielectric permitiv- 
ity in C 

In this section we establish Theorem l2.2l Recall that any element u oi H^^^{Y \R) / C can be written 
as 

u = 7'iu + 7'2U+ ^ Vx^u + d. (6.1) 

-i<A„<i 

where d is chosen such that fyy^jiudy = 0. Here Vi are the orthgonal projections onto Wi, i ^ 1,2 
and V\^^ are the orthogonal projections onto the finite dimensional eigenspaces associated with the 
eigenvalues A„ in (—1/2, 1/2). The orthogonal decomposition (|6.1I) is used in the proof of Theorem 
given below. 



Proof. For u, v in iJpg^(F \R)/C we apply (|6.ip to see that 

2 

= (7'1U,W) + Z(P2M,W)+ X! (l + (-2-l)(^-^"))(^A„?^,v) 

-i<A„<i 

= (T,u,tO. (6.2) 
From (|6.2p we conclude that 

=7'i+z7'2+ ^ (l + (z-l)(i-A„))7'A„. (6.3) 

-i<A„<i 

It is evident from (|6.3p that for z ^ (A; + l/2)/(Ai — 1/2) that Tz is a bounded one to one and onto 
map in H^^^iY \ R)/C The formula for T^^ is given by 

^V,+z-'r2+ (l + (^-l)4-^™))"'^A„. (6.4) 

-i<A„<i 



Taking conjugates on both sides of (16.21) gives Bz{u,v) — {v,Tzu) and choosing u = T^^q for q in 
Hl^^{Y \ R)/C delivers the identity 



Bz{Tz7'q,v)^{v,q). (6.5) 

To complete the proof consider any linear functional F in [Hp^j.{Y \ R)/C\* with F{v) = for 
V = const. Applying the Reisz representation theorem shows that there exists a unique solution u 
inHl,,{Y\R)/C of 



Bz{u, v) = F{v), for all v in H^^^{Y \ R)/C, (6.6) 

□ 

and Theorem 12.21 is proved. 
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7 Convergence of the Power Series 

In this section we show that the series (|4.ip . (|4.4p identified in Theorems 14.11 and 14.21 are convergent. 
Here the convergence radius R depends upon the branch of the quasistatic dispersion relation = 
^0 (t, k). The methodology applied here uses generating functions and follows the approach presented 
in [24]. 

Step I. A priori estimates. 

We establish a priori bounds for the solutions of p.2p and (|3.15l) . Before proceeding recall the 
Friderich's inequality and trace estimate satisfied by elements V of Hp^j.{Y \R)/C given by 

My\r= (^J^^^mi^ + m'dy^ <n\m (7.1) 

and 

IIV'II <AU\\ (7.2) 

where Q, and A are positive constants depending on R and 11-011 is the norm on Hp^^{Y\R)/<C defined 
by (1221). 

The solutions tp of (|3.2p solve boundary value problems of the following generic form. Given 
ip e H^^^IJ \R),Ge L^{R) and F e [L^iR)]^ the function is the H^{R) solution of 



( V0 • VzJ - ^o^v) = ^ Gv + Jj^F-\/v yv e (R) 



(7.3) 



Decompose tfj according to tp — ipa + ipi + 'tp2 such that 



/«(VVo-VTJ-eo^oi^) = /^JG+(l+eo)^i + (l + eo)^2]^J yveH^iR) 
i'oldR^ = 0, 



j /^(V0i • + 01 w) = Vv e (R) 

[i'lldR- = 0|9fl+, 

and 

f /fl ( Vi^2 • VzJ + V^2i^) = F • VzJ Vt; e 1 (i?) 

1 02|9_R- = 0. 



(7.5) 



(7.6) 



Let J = G + (1 + ^o)V'i + (1 + ^0)02 denote the right hand side for (|7.4p . Expanding ipQ and J with 
respect to the complete orthonormal set of Dirichlet eigenfunctions {4>j}°^i, gives 



\MmiR.)^T.(^ 



(Co - l^j) 



<ClJJ\\h^n)^ (7.7) 



where 



(1 + ^,)V2 



Gi,„ = max <^ maxj < — — } } . (7. 



Collecting results gives 



Wi^iWmiR) = \\i^2\\m/2^aB) < 

||02||ffi(ii) < ||^'||L2(i^) 

WMhhr) < Ci,„\\G + {1 + + (l + OMlhr) (7.9) 
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and 

miHHR) < i3„(||G|U2(^) + ||F|U.(^) + ll^ll), (7.10) 

where 

B„ = mai5„e/„{Ci,„, Ci,„(l + fo), ^(Ci,„(l + Co) + !)}• (7.11) 
Now we estimate the field V'* • The explicit expression p.6p for ^» gives 

and 

\\tp*\\H^R) < Ln- (7.13) 

where L„ is the maximum of the right hand side of (|7.12p for (ti '*) G In ■ 

Next we provide an upper bound for the solutions of ip'^ of (j3.15l) . We may continuously extend 
elements v in iJ^g,. (F \ i?)/C onto R as Hp^^iY) functions such that the extension u*^^* satisfies 

K'-^'Wy < C\\v\\. (7.14) 

In what follows we continue to denote these extensions by v and (|3.15l) takes the equivalent form 
p.Sp . Application of Holders inequality to the right hand side of ()3.8p gives 

T^\BMm,v)\ = C {\\F^\\L-t^Y\B) + \\Gi\\l-[y\b) + \\F2\\l-[R) + \\G2\\l-[R)) MV (7.15) 

It is evident that Bziipm^'^'m) ~ (^zV'miV'm) ^^'^ applying the explicit representation for Tz 
delivers the estimate 

< CxGz {\\Fi\\l2^y\r) + \\Gi\\L2iY\R) + \\F2\\mR) + I|G2||l2(^)) ||z;||, (7.16) 

with 

= (min{l, |z|, _^/2i,nf ^^^^{|1 + (z - l)(l/2 - Xn)\}})~'. (7.17) 
For z = ep(Co) we maximize (|7.17p for G -^n and choose v = ifj'^-^ in (|7.16p to obtain 

AWJ] < G„ {\\f,\\l2^y\r} + \\Gi\\mY\R) + mWmR) + \\G2\\mR)) ■ (7.i8) 

Applying Freidrich's inequality we arrive at the desired upper bound given by 

'^^IIV'mlll'\i?. < -E-n {\\Fi\\l^{Y\R^) + 1 1 Gi 1 1 ^2 + | I-F2 | U2(j^,) + 1 1 G2 1 1 L2 (fl.) ) , (7.19) 

where £"„ denotes a generic constant depending only on /„. 

The a priori estimate for ip follows from the representation formula (|3.18p and 

ll^l|y\fl < (7.20) 

where the constant _ff„ depends only on /„. 
Step II. System of inequalities. 
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We apply the a priori estimates developed in Step I to the system of equations (|3.2I) , (|3.15p , and 
solvability conditions (I3.19|) derived in section [31 From (13.81) and (I7.19|) . it follows that 



-^\Wm\\Y\R< 



m-2 



En[T\Yl \il\U^-l\\Y\B+Y. \il\U^-l-l\\Y\R+Yl M\'P^-2-i\\y\R 
1=1 1=1 1=1 

m — 2 m — 3 

+ Mm-l\\Y\R + \\lpm-2\\Y\R) + 1^' 1 1 1 V'm-2-/ 1 li/i (i?^) + 16 1 1 1 V'm-3-/ 1 1 ifi («) 

m — 4 

+ Y \^l\\\'^m-4-l\\H^iR) + \\lpm-2\\m{R) + I IV'm-3 I («) + 1 1 V'm-4 1 1 ffi (i?) 
1=1 

m-2 I 

-n\\m{R)) 

1=0 n=0 
m-2 

+ E MU^-2-l\\Y\R+Um-2-l\\m{R))]- (7.21) 



/=0 

From and (fTTU)) . it follows that 

W-lpmWmiR.) < 



-1 



Bn[Y ilW'^rn-lWHHR) + ^ ^ 1 1 V'm- 1-i 1 1 ^ (fi) 
1=1 1=0 
m-2 

+ E C'IIV'm-2-i||Hi(«) + IIV'm-llkM^) + IIV'm-2||i/i(fl.) 
1=0 

m— 1 I m—1 

+ E E6"-'^«ll^''-«II^Mfl) + E ^IIV-^-'IUm^) + (7-22) 

i=l n=Q 1 = 1 

The solvability constant Q in p.l9p is bounded away from zero and oo for G In Put m — 2 i— m 
in p. 191) to obtain the inequality for The inequality is in terms of a constant C„ depending 

only on /„ and is given by 

\u\ < 

m—1 771—1 
Cn[r'{ E M\i'm+l-l\\Y\R + ^ M\^n^-l\\Y\R + 2||C+lllnfl + \My\r) 

m — 1 m — 2 

+ E l^'IIIV'm-l-zllffi(fl) + E l^nilV'm-2-z||//i(^?) 
1=0 1=0 

m — 1 I m—1 

+ E E \^-m~l\\U{\\'4>l-n\\miR) + ||Vi-«||Y\fl) + 2 ^ |6 | (| I^Am-i | + 1 1 V'm-i 1 1 ) 

1=1 n=0 1=1 

HUmWmiR) + \\^m\\Y\R) + ( 1 1 V^m 1 1 i (i?,) + llV-mllnfl)] (7-23) 

The decompositions and bounds ([7^31) . (fT?^ give 

ll'/'mllffM^) - and (7.24) 

Um\\Y\R < U'm\\Y\R + H,,\U-ll (7.25) 

Here the upper bounds ([7^ . ((7?^ ((7?^ . (17^ . and hold for all values of r for which 

the branch of the quasistatic dispersion relation lies in interval /„. 
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Step III. Majorizing sequence and analyticity of generating functions. 
We simplify the exposition by introducing 



Pjn 




\y\r. 


Pm — 


||T"Vm 


\h^(r) 


Pm — 






Pm 




\y\r 


— 







(7.26) 
(7.27) 
(7.28) 
(7.29) 
(7.30) 

and show that the corresponding series J2p"^Pm: Y^P^Pm, ^P^^Pm, J2p"^Pm: ^^'^ J2p"^^rn con- 
verge. This is sufficient to estabhsh summabihty and convergence for the series (|4.4p and (j4.ip . 
Writing (fT!2T|) . fTT^ (fTMl) . ([7^241) . and (f7!25l) in terms of the new notation gives the following 
system of inequalities. For to > 2 we have: 

pL < 

771 — 2 m — 2 m — 2 

En [ Sip.m_i + T ^ SW^__^_i + ^ Sfp„,_2_i 

1=1 1=0 1=0 

m — 2 m — 3 m— 4 

+ 'r^Pm-2 + X! ^lP^n-2-l +"^51 ^lP"i-3-l + '^^ ^ SiPjn-4-l 

1=0 1=0 



1=0 



m-2 I 



'Pm-2 + TPm-3 + + X! X! ^rn,-2~lSniPl-n + Pl-n) 

1=0 n=0 

m-2 



«'(Pm-2-( +J5m-2-i)], 



1=0 



and for to > 1, 

Pm < 



Bn [ Y, ^iPrn-l + T Y^ Sip„i-i-l + ^ SlPm-2-l + Tpm-1 + T^Pm-2 
1=1 1=0 1=0 

m — 1 I m— 1 

+ ^ ^ Sm-lSnPl-n + ^ S;P,„_; + p„] . 



1=1 n=0 1=1 

The bounds ((7?^ . (TOi)) . and ((7?^ yield the following bounds for to > 1: 

m — 1 m — 1 

^2 



m — 1 rn — 2 m — Z I 

+T Y SlP,n-l-l + X! *'P™-2-i + X! X! *™-iSn(P/_„ + P/-n) 
1=0 1=0 1=0 n=0 

m 

+ 2 5Z^'(^™-i +P^-l) + + (2 + T^)Pm]: 



1=1 



Pm — Pm ~^ LnSrm 
Pm ^ Pm+ 



(7.31) 



(7.32) 



(7.33) 

(7.34) 
(7.35) 
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The inequalities presented above hold for all r for which belongs to /„. In what follows it is 

convenient to increase the upper bounds by choosing the maximum value of r for in 7„ and 
absorbing this value into the constants Bn,Cn,En Additionally any constant factors other than 

unity that multiply terms in these upper bounds are absorbed into these constants. With these 
choices the system is written: 

P'm < (7.36) 

m-2 



m-2 

SlPm-l- 


m-2 


I 


1=0 


1=0 




m-2 


m-3 


m-4 


Y SlPm-2-l 
1=0 


: + X] ^lPm-3-l 4 

;=o 


- X ^lPm-4-. 
1=0 



m-2 



+Pm-2 + Pm-S + Pm-4 + X/ X/ ^rn-2-lSn{Pl-n + Pl-n) 

1=0 n=0 

m-2 

+ X ^l{Pm-2-l +Pm-2-l)\ , 



Pm < (7.37) 

m — 1 m— 1 m— 2 

Bn[ ^ SlPm-l + ^ SiPm-1-l + ^ SlPm-2-l +Pm-1 +l>m-2 

/=o /=o 

m — 1 / m — 1 

+ X X ^rn-lSnPl-n + X ^iPm-l +Pm] > 
1=1 n=0 1=1 



Sm < (7.38) 

m — 1 m — 1 

Cn[Pm+l + X ^iPm+l-l + X ^iPm-l 
1=1 1=0 
m — 1 m— 2 m — 3 / 

+ X ^iPm-l-l + X ^lPm-2-l + X X «™-i«"(P(-n + Pl-n) 
1=0 1=0 1=0 n=0 



X (Pm-l + Pm-l) +Pm+Pr, 



Pm < Ci{pm + Sm), (7.39) 
Pm<C5ip'm+Sm-l), (7.40) 



where C4, C5 are positive constants depending only on /„ and m > 1. 
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We now introduce a majorizing sequence {{am, ^m, Cm.dmi em)}m^o which 

oo oo 

^ a„p" = pa{p) +p„ >Y^ p„p" 



n=0 n=0 
oo oo 

Y.bnP"= b{p) >Y.PnP" 
n—1 n—l 

oo oo 



n=0 n=0 

oo oo 

oo oo 

^e„p"= pe(p) >5^P>", (7.41) 



n=2 n=2 

to show that the generating functions a(p), . . . ,e(p) are analytic in a neighborhood of the origin 
p = 0. 

The majorizing sequence is chosen so that the system of inequahties (|7.36p . (|7.37p . (|7.38l) . (|7.39l) . 
(|7.40p hold with equality for = &m = Pm, Cm = Sm, Pm = dm, and = p^^. Indeed, for 
this choice one observes that p„j < am, Pm ^ bm , Sm < Cm, Pm < dm, and < e„. Enforcing 
equality in (|7.36l) . (|7.37l) . (|7.38l) . (|7.39l) . (|7.40|) . multiplying each by the appropriate power of p 
and summation delivers the equivalent system for the generating functions a{p), b{p), c{p), d{p), e{p) 
given by: 

A^{a, b, c, d, e, p) = 0, i = 1, . . . 5, (7.42) 

where: 

Ai{a,b,c,d,e,p) = p^ - a + Csie + {c - so)), (7.43) 

A2{a,b,c,d,e,p)= (7.44) 
-6 + S„[(c- so){d-po) + {cd + d){p + p^) + c^d- Sppo - psoSi(Po + Pi) 
-soPo(c - So) - SQcrf + so(p(soPi + siPo) + soPo) + (c - so)(c' - Po) + pa] 

A^{a,b,c,d,e,p) ^ (7.45) 
-(c - So) + Cn [(a — Pi)(c — So) + e + pac + p(a + cd) 
+p^cd + c^(pa +po) - soPqC - so((pa +Po)c - soPo) + (c - so)(d - pa) + h 
-pipa+Po)csi - p^S2(pa + pQ)c + c^d- sqPqc ~ so(cd- sopo) - pcd(si + ps2)] 

A4a, b, c, d, e, p) = (d - Po) + C^ib + (c - so)). (7.46) 

A5(a,&,c,d,e,p) = (7.47) 
-e + i?„[(a - Pi){c - So) + pac + p{pa + pQ)c + pa + p{pa + pQ)) 
+(c + l)d{p + p^ + p^) + pc{c + l)(pa + Po + d)] 

One can check Ai, A2, A3, A4 and As vanish for a — Pi — \\^i\\y\r, b — 0, c = sq = d = pq ^ 
\\4'o\\ {R) , e = 0, and p = 0. Moreover, the determinant of the Jacobian matrix with respect to the 
first five variables at this point, i.e., (a, 6, c, d, e, p) — (p^, 0, so,po, 0, 0), is 

det —--r — :r~\ Pi>0>so,po,0,0 = 1. 7.48 

o(a, 0, c, a, e) 
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Thus the impHcit function theorem for functions of several complex variables guarantees that the 
generating functions a, h, c, d, e are analytic within a neighborhood of p = and we conclude that the 
series on the right hand side of the inequalities (|7.41l) are convergent. Here pp = llV^ollrVfl = 
and the initial data poj Pi £^re bounded functions of on /„. Choosing the maximum values of poj 
Pi^ ^0 1 Co varies over /„ delivers a majorizing sequence with radius of convergence R depending 
only on J„. It now follows that the series (|4.1[) converges and that ()4.4p is summable for < p < i?. 

8 Power Series Solution of Cell Problem 

In this section we show that the functions defined by the power series solve l|2.4p . For in /„ recall 
the series ([2:6| 



u — 

m=0 
oo 



V'^Um (8.1) 

^ V'^U (8.2) 



m=0 



with Urn = i^'UQipm and ij — pr converges for < p < i?. For v G -ffper(^) < p < R. We 
define 

a''(w) := / T^{(, ~ er^){'Vy + ir]k)u ■ (Wy + ir]k)v + / T^^{Vy + ir]k)u ■ {Wy + if]k)v 
Jh c Jp 

7f{^- er^){Wy+irjk)u- {Vy+iT]k)v- / i]^ £,{£,- er^)uv. (8.3) 

Jy 

First observe that aP{v) has a convergent power series in rj that is obtained by inserting u and ^ 
into the above expression and expanding in powers of rj . The coefficients of this expansion are 
exactly the left-hand side of equation (12.81) . Since ijjm satisfies (|2.8I) . all coefficients of the power 
series expansion of a'^{v) are zero and we conclude that a'^{v) = for all v G Hp^j.{Y), < p < R. 
Thus we conclude that the pair (m, C) satisfies (|2.4p . 

Next we record the following property of the sequence {£,m} that follows directly from the vari- 
ational formulation of the problem and the power series solution. 

Theorem 8.1. The functions C„i cire real for all m. 

Proof. Setting v = u in (|2.4p delivers a quadratic equation for ^. Since the discriminant greater 
than zero, we conclude that ^ is real and it follows that ^„i, m — 1, • . ., are real. □ 

With these results in hand, Theorems 14.11 and 14.21 now follow from the convergence established in 
section [7] together with Theorem 13. II 



9 Acknowledgments 

This research is supported by NSF grant DMS-0807265 and AFOSR grant FA9550-05-0008. 



References 

[1] Alu, A., Engheta, N: Dynamical theory of artificial optical magnetism produced by rings of 
plasmonic nanoparticles. Phys. Rev. B Vol 78, 085112 (2008) 

[2] Adams, D.: Sobolev Spaces. Elsevier, Waltham(2003) 



25 



Ammari, H., Kang, H., Lee, H.: Asymptotic analysis of high-contrast phononic crystals and a 
criterion for the band-gap opening. Arch. Rational Mech. Anal. Vol 193, 679-714 (2009) 

Ammari, H.,Kang, H., Soussi, S., Zribi, H.: Layer potential techniques in spectral analysis. Part 
IL Sensitivity analysis of spectral properties of high contrast band-gap materials. SIAM Multi. 
Model. Simu. Vol 5, 646-663 (2006) 

Bergman,D.J.: The dielectric constant of a simple cubic array of identical spheres. J. Phys. C 
Vol 12, 4947-4960 (1979) 

Bohrcn, C.F., Huffman, D.H.: Absorption and Scattering of Light by Small Particles. Wiley, 

Wcinheim (2004) 

Bouchitte, G., Bourel, C. : Homogcnization of finite metallic fibers and SD-effective permittivity 
tensor. Proc. SPIE, Vol 7029, 702914 (2008); doi:10.1117/12. 794935 

Bouchitte, G., Schweizer, B.: Homogenization of Maxwell's equations in a split ring geometry. 
SIAM Multi. Model. Simu. Vol 8(3), 717-750 (2010) 

Bouchitte, G. , Felbacq, D.: Negative refraction in periodic and random photonic crystals. New 
J. Phys. Vol 7, 159(2005) 

Bouchitte, G. , Felbacq, D.: Homogenization near resonances and artificial magnetism from 
dielectrics . C. R. Acad. Sci. Paris 1(339), 377-382 (2004) 

Chen, Y., Lipton, R.: Tunable double negative band structure from non- magnetic coated rods. 
Now .J. Phys. Vol 12, 083010 (2010) 

Cherdantsev, M.: Spectral convergence for high-contrast elliptic periodic problems with a defect 
via homogenization. Mathematika Vol 55, 29-57 (2009) 

Chern, R. L., Felbacq, D.: Artificial magnetism and anticrossing interaction in photonic crystals 
and split-ring structures . Phys. Rev. B Vol 79(7), 075118 (2009) 

Colton, D. , Kress, R.: Inverse Acoustic and Electromagnetic Scattering Theory. Springer, New 
York (1998) 

Dolling, G., Enrich, C., Wegener, M., Soukoulis, C. M., Linden, S.: Low-loss negative-index 
metamaterial at telecommunication wavelengths. Opt. Lett. Vol 31, 1800-1802(2006) 

Felbacq, D., Bouchitte, G.: Homogenization of wire mesh photonic crystals embdedded in a 
medium with a negative permeability. Phys. Rev. Lett. Vol 94(18), 183902(2005) 

Felbacq, D., Guizal, B., Bouchitte, G., Bourel, G.: Resonant homogenization of a dielectric 
metamaterial. Microwave Opt. Tech. Lett. Vol 51(11) 2695-2701(2009) 

Farhat, M., Guenneau, Enoch, S., Movchan, A. B.: Negative refraction, surface modes, and 
superlensing effect via homogenization near resonances for a fi.nite array of split-ring resonators. 
Phys. Rev. E Vol 80(4), 046309 (2009) 

Figotin, A. , Kuchment, P.: Band-Gap Structure of the Spectrum of Periodic Dielectric and 
Acoustic Media. I. Scalar model. SIAM J. Appl. Math. Vol 56(1), 68-88 (1996) 

Figotin, A. , Kuchment, P.: Band-Gap Structure of the Spectrum of Periodic Dielectric and 
Acoustic Media. II. 2D Photonic Crystals. SIAM J. Appl. Math. Vol 56(6), 1561-1620 (1996) 

Figotin, A. , Kuchment, P.: Spectral properties of classical waves in high contrastperiodic 
periodic media. SIAM J. Appl. Math. Vol 58(2), 683-702(1998) 



26 



[22] Folland , G. : Introduction to Partial Differential Equations. Princeton University Press, Prince- 
ton (1995) 

[23] Fortes, S. P., Lipton, R. P., Shipman, S. P.: Sub-wavclcngth plasmonic crystals: dispersion 
relations and effective properties. Proc. R. Soc. Lond. Ser. A Mat. doi: 10. 1098/rspa. 2009. 0542 
(2009) 

[24] Fortes, S. P., Lipton, R. P., Shipman, S. P.: Convergent power series for fields in positive or 
negative high-contrast periodic media. Comm. Partial Differential Equations Vol 36(6), 1016- 

1043(2011) 

[25] Hempel, R., Lienau, K.: Spectral properties of periodic media in the large coupling limit. 
Comm. Partial Differential Equations Vol 25, 1445-1470 (2000) 

[26] Huangfu, J., Ran, L., Chen, H., Zhang, X., Chen, K., Grzegorczyk, T. M., Kong, J., A.: 
Experimental confirmation of negative refractive index of a metamaterial composed of f2-like 
metaUic patterns. Appl. Phys. Lett. Vol 84, 1537(2004) 

[27] Huang, K., C, PovincUi, M., L., Joannopoulos, J. D.: Negative effective permeability in polari- 
tonic photonic crystals. Appl. Phys. Lett. Vol 85(4), 543 (2004) 

[28] Jelinek, L. , Marques, R.: Artificial magnetism and left-handed media from dielectric rings and 
rods . J. Phys.:Condens. Matter Vol 22(2), 025902(2010) 

[29] Kamotski, V., Matthies, M., Smyshlyaev, V.: Exponential homogenization of linear second 
order elliptic pdes with periodic coefficients. SIAM J. Math. Anal. Vol 38(5), 1565-1587 (2007) 

[30] Kohn, R., Shipman, S.: Magnetism and homogenization of micro- resonators. SIAM Multiscale 
Model Simul. Vol 7(1), 62-92 (2008) 

[31] Lepctit, T. ,Akniansoy, E.: Magnetism in high-contrast dielectric photonic crystals. Microwave 
Opt. Tech. Lett. Vol 50 (4), 909-911(2008) 

[32] Milton, G. W.: Realizability of metamaterials with prescribed electric permittivity and magnetic 
permeability tensors. New J. Phys. Vol 12(3), 033035(2010) 

[33] Milton, G.W. :The Theory of Composites. Cambridge University Press, Cambridgc(2002) 

[34] O'Brien, S., Pendry, J. B.: Photonic band-gap effects and magnetic activity in dielectric com- 
posites. J. Phys.: Condens. Matter Vol 14, 4035-4044(2002) 

[35] Pendry, J., Holden, A., Robbins, D., Stewart, W.: Magnetisim from conductors and enhanced 
nonlinear phenomena. IEEE Trans. Microwave Theory Tech. Vol. 47(11), 2075-2084 (1999) 

[36] Pendry, J., Holden, A., Robbins, D., Stewart, W.: Low frequency plasmons in thin-wire struc- 
tures. J. Phys.: Condens. Matter Vol 10, 4785-4809 (1998) 

[37] Peng, L., Ran, L., Chen, H. , Zhang, H., Kong, L. A., Grzegorczyk, T. M.: Experimental 
observation of left-handed behavior in an array of standard dielectric resonators. Phys. Rev. 
Lett. Vol 98, 157403 (2007) 

[38] Shvets, G., Urzhumov, Y.: Engineering the electromagnetic properties of periodic nanostruc- 
tures using electrostatic resonances. Phys. Rev. Lett. Vol 93(24), 243902-1-4 (2004) 

[39] Shipman, S.: Power series for waves in micro-resonator arrays. In: Proceedings of the 13th 
International Conference on Mathematical Methods in Electrodynamic Theory, MMET 10, 
Kyiv, Ukraine September 6-8, 2010, IEEE (2010) 



27 



[40] Veselago, V., G.: The electrodynamics of substances with simultaneously negative values of e 
and fi. Sov. Phys. Usp. Vol 10, 509 (1968) 

[41] Vynck, K., Felbacq, D., Centeno, E., Cabuz, A., I., Cassagne, D., Guizal, B.: All-dielectric 
rod-type metamaterials at optical frequencies. Phys. Rev. Lett. Vol 102, 133901 (2009) 

[42] Service, R. F.: Next Wave of metamaterials hopes to fuel the revolution. Science Vol 327 
(5962), 138-139(2010) 

[43] Shalaev, V.: Optical negative-index metamaterials. Nature Photonics Vol 1, 41-48 (2007) 

[44] Shalaev, V. M., Cai, W., Chettiar,U. K., Yuan, H. K., Sarychev, A. K., Drachev, V. P., Kild- 
ishev, A. V.: Negative index of refraction in optical metamaterials. Opt. Lett. Vol 30(24), 
3356-3358(2005) 

[45] Smith, D., Padilla, W., Vier, D., Nemat-Nasser, S. Schultz, S.: Composite medium with simul- 
taneously negative permeability and permittivity. Phys. Rev. Lett. Vol 84 (18) 4184 4187(2000) 

[46] Wheeler, M. S., Aitchison, J. S., Mojahedi, M. : Coated non-magnetic spheres with a negative 
index of refraction at infrared frequencies . Phys. Rev. B Vol 73, 045105(2006) 

[47] Yannopapas, V.: Negative refractive index in the near-UV from Au-coated CuCl nanoparticle 
superlattices. Phys. Stat. Sol. (RRL) 1(5), 208-210(2007) 

[48] Yannopapas, V.: Artificial magnetism and negative refractive index in three-dimensional meta- 
materials of spherical particles at near-infrared and visible frequencies. Appl. Phys. A Vol 87 

(2), 259-264(2007) 

[49] Zhang, P., Potet, S., Carbonell, J., Lheurette, E., Vanbesien, O., Zhao, X., Lippens, D.: 
Negative-zero-positive refractive index in a prism-like omega-type metamaterial . IEEE Trans. 
Microw. Theory Tech. Vol 56, 2566(2008) 

[50] Zhang, S., Fan, W., Minhas, B. K., Frauenglass, A., Malloy, K. J., Brueck, S. R. J.: Midinfrared 
resonant magnetic nanostructures exhibiting a negative permeability. Phys. Rev. Lett. Vol 
94,037402(2005) 

[51] Zhikov, V. v.: Gaps in the spectrum of some elliptic operators in divergent form with pe- 
riodic coefficients. Algebra i Analiz Vol 16 34-58 (2004). English translation St. Petersburg 
Mathematical Journal Vol 16(5) 773-790 (2005) 

[52] Zhou, X., Zhao, X. P.: Resonant condition of unitary dendritic structure with overlapping 
negative permittivity and permeability. Appl. Phys. Lett. Vol 91, 181908(2007) 



28 



